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The Extended Theory of Finite Fermi Systems (ETFFS) describes nuclear excitations considering 
phonons and pairing degrees of freedom, using the effective Landau-Migdal interaction and nuclear 
mean fields obtained from experimental data. Here we employ the nuclear mean field derived from 
Skyrme interactions and the corresponding particle hole interaction. This allows to extend the range 
of applicability of the ETFFS to experimentally not yet investigated short-lived isotopes. We find 
that Skyrme interactions which reproduce at the mean field level both ground state properties and 
nuclear excitations are able to describe the spreading widths of the giant resonances in the new 
approach, but produce shifts of the centroid energies. A renormalization of the Skyrme interactions 
is required for approaches going beyond the mean field level. 



I. INTRODUCTION 

The renewed interest in nuclear structure physics is 
due to its appUcation to astrophysics. Theories for the 
element synthesis in stars need both bulk properties of 
unstable isotopes like masses, and more specific informa- 
tion, such as level densities and neutrino-nucleus cross 
sections. Eventually, the compression modulus and the 
symmetry energy of neutron-rich nuclei will be deter- 
mined. For a recent review, see e.g. Ref. 

Most calculations in this field are restricted to the 
mean field approximation which describes excited states 
by the one-particle one-hole random phase approxima- 
tion (Iplh RPA). In this approximation one is able to 
describe two integral characteristics of giant resonances: 
the mean energy and the integrated strength. The width 
of the giant resonances is underestimated in mean field 
approaches, however, even if one includes the single par- 
ticle continuum. 

In unstable nuclei, the analysis of, for example, the 
breathing mode and the giant electric dipole resonance 
requires realistic strength distributions because these res- 
onances are expected to be broad and various multipole 
resonances overlap one another. Theories beyond the 
mean field level are called for. 

The Extended Theory of Finite Fermi Systems 
(ETFFS) is a linear response theory which incorporates 
phonons and pairing degrees of freedom. One essentially 
replaces the particle-hole propagator, which leads to the 
random phase approximation, by a more general prop- 
agator and changes, correspondingly, the renormalized 
interaction, the effective operators and the ground-state 
correlations. Details can be found in Refs. [1, 0, H, Q, 
where the authors worked out the coupling of phonons 
to the single particle propagators and to the interac- 
tion in a consistent way. The present calculations are 
based on the Quasi particle Time Blocking Approxima- 
tion (QTBA) that also considers pairing correlations 7]. 
Examples and references to previous calculations and mi- 
croscopic analyses of experimental data can be found in 



Another model, the QRPA plus phonon coupling 
model, has been developed by Bortignon and coworkers 

Most excited states predicted by the conventional 
random-phase approximation depend sensitively on the 
level density of the single-particle spectrum, as is well 
known from the schematic model by Brown and Bolsterli 
[ni- As the TFFS, with and without pairing, uses the 
experimental single particle spectrum as far as it is avail- 
able, e.g. it generally comes closer to the ex- 
perimental data than self-consistent approaches do. For 
that reason an extended version of the TFFS in which 
the effects of the phonons are included is the appropriate 
tool for systematic investigations of both collective and 
non collective states [ij]. Its applicability, however, is 
limited by the need for experimental information on the 
level density in the neighboring odd mass nuclei. 

At present, these properties of most short-lived neu- 
tron rich nuclei are not known experimentally. Therefore 
one is forced to use nuclear mean fields extrapolated from 
the known stable nuclei. One of the most successful ap- 
proaches starts with Skyrme forces, a class of effective 
interactions defined at the mean field level. When going 
beyond the mean field level, one has to expect a renor- 
malization of the interaction strengths p^ . 

In the present communication, we investigate quanti- 
tatively to which extent a renormalization of the Skyrme 
parameterization may be required by employing a self- 
consistent mean field and the corresponding particle hole 
interaction derived from existing Skyrme forces in the 
framework of the Extended Theory of Finite Fermi sys- 
tems. 

There are various Skyrme parameterizations that re- 
produce the ground state properties of the known iso- 
topes with similar accuracy but result in different mean 
fields, as reviewed recently in Refs. .161 [l7t. 

The nuclear matter properties of the Skyrme param- 
eterizations, such as the effective mass, the compression 
modulus, and the symmetry energy can be used to group 
the effective interactions into some broad classes. The 
versions initiated by Vautherin and Brink were fitted to 
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ground state properties only. Those interactions are char- 
acterized by a large incompressibility (see SIII in Table 
H]). Later versions of the Skyrme interactions modified the 

TABLE I: Nuclear matter properties of the Skyrme param- 
eterizations used in the present calculation; for comparison 
also SIII is given 
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density dependence of the interaction in order to allow a 
smaller incompressibility. Those interactions reproduce 
the mean excitation energies of giant resonances reason- 
ably well and are therefore a natural starting point for 
our investigation. 

The interaction SLy4 ^19*1 has been used in(Q)RPA 
[20l . [21I [2^ and in some first attempts to include explic- 
itly long range quadrupole correlations {alias quadrupole 
phonons) into the nuclear ground state [23j . An inves- 
tigation similar to the present one, but based on the 
relativistic mean field approach, has been published re- 
cently (23.[25l|. Their results, calculated with and without 
phonons differ very little from each other, in contrast to 
our results. The reason is the so-called subtraction pro- 
cedure of Ref. which these authors apply. This is an 
ad hoc condition that is difficult to justify, as one needs 
a consistent description of both even even and even odd 
nuclei. 



II. METHOD 

We calculate the mean field within the HFB approach 
using the code HFBRAD [2^] and the excited states with 
the the QTBA version 0] of the ETFFS. Here, as men- 
tioned before, the authors developed a renormalization 
procedure that corrects for possible double counting of 
the phonon contributions that are already included in the 
Skyrme parameterization. We leave this correction out. 
This is also in line with the investigations by Bertsch et 
al. [2^, who argue that an effective interaction has to be 
refitted if one includes additional correlations. As usual, 
the residual interaction for (Q)RPA and the present the- 
ory is derived as the second derivative of the Skyrme 
energy functional [13, ■ In the present calculation we 
make some simplifications. For the non magic nuclei we 
solved the HFB equations and extracted the quasi parti- 
cle energies, wave functions and the occupation numbers. 
In the second derivative, the ph-interaction, we dropped 
the spin-orbit term and we approximated the velocity- 
dependent terms by its Landau-Migdal limit [13, |28| . 
There are two kinds of velocity-dependent terms: the 
first is oc k^(5(r — r') and the second one is oc k^^(5(r — r')k 



(p-wave interaction in momentum space). The value of 
the first term averaged over the angular variables gives 
^kp6{r — r'), while that of the second one vanishes. Such 
an approximation, of course, violates self-consistency and 
one has to modify the parameters of the residual interac- 
tion to obtain the spurious ccnter-of-mass state at zero 
energy. For this reason we changed the parameters pro- 
portional to tikp6{r — r') up to 20%. 

In principle, the ETFFS treats the single-particle con- 
tinuum on the RPA level for magic nuclei exactly and 
includes the ground state correlations caused by CC 
(GSCp/ion) [1] explicitly. However, because of techni- 
cal difficulties caused by the pairing correlations, we did 
not consider these effects in the present calculation. The 
ground state correlations are therefore restricted to the 
conventional RPA correlations. We discretized the con- 
tinuum with a quasi-particle energy cutoff of 100 MeV. 
We did not find any noticeable differences in the results 
with a larger basis. We used 14-16 low-lying phonons 
with angular momentum L = 2 — 6 and natural parity. 



III. RESULTS 

A. Collective isovector electric dipole states 

Here we present our results for the giant dipole ex- 
citation for some of the tin isotopes from ^'^°Sn up to 
^^^Sn, where we used the SLy4 [l9| parameterization of 
the Skyrme force. The SLy4 [19| parameterization is not 
only adjusted to ground state properties, but it also re- 
produces within the RPA the GDR and isoscalar collec- 
tive states reasonably well. For comparison we used for 
^^^Sn two further parameterizations, namely BSk5 [1^ 
and SKM* [s^l . We compare our theoretical results with 
the data of ^^°Sn and ^■^^Sn. In general, the experimental 
GDR data are analysed by fitting to a Lorentzian. For 
that reason we have used the same procedure to deter- 
mine the calculated integral characteristics [^, [3l| . The 
mean energies Eq and the resonance widths F have been 
extracted from the calculated strength distribution under 
the condition that the first three energy-weighted mo- 
ments of the Lorentzian and the theoretical distribution 
of the total photo absorption cross section should coin- 
cide. This method works well in the case of the ETFFS 
results. 

In Fig. [T]the photo absorption cross sections for ^^°Sn 
and ^'^■^Sn are compared with the experimental ones. The 
full line represents the result of the present approach. 
The dotted line is a Lorentz parametrization of the ex- 
perimental data. In both cases the shape of the calcu- 
lated cross sections agrees very well with the experimen- 
tal ones. The energies are in both cases, however, too 
low. Quantitatively, we give the same results in Table HIl 
and compare it with the RPA results. In Figl2]the exper- 
imental cross section of ^'^^Sn (dotted line) is compared 
with the results of the RPA and the present theory where 
we used three different Skyrme parameterizations. The 
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FIG. 1: El photo absorption cross section for Sn and Sn 
calculated within the present theory (full). The dotted lines 
indicate the data for ^^°Sn [H] and "^Sn The smearing 
parameter is A = 0.2 MeV. The SLy4 Skyrme parametriza- 
tion was used. 



theoretical results look similar but deviate in detail. In 
all three cases the energies calculated within the ETFFS 
approach are too low compared to the data (s^l , but the 
strength distribution is strongly improved compared to 
the RPA calculation. In Table |TT] we give a more detail 
comparison. The BSk5 and SLy4 parameterizations give 
very similar results; the widths are roughly 50% larger 
compared to the one calculated in RPA and are in good 
agreement with the data. The energies of the extended 
theory are more than 1.5 MeV lower than the RPA re- 
sults which are already 2 MeV too low compared to the 
experiment. 

Both effects are connected with the inclusion of the 
phonons. (I) The single particle strength is fragmented 
due to the coupling to the phonons and this effect gives 
rise to the spreading width which strongly improves the 
agreement with the experiments. (II) The phonons, how- 
ever, also compress the single particle spectrum and for 
that reason the energies calculated in ETFFS are lower 
by about 1.5 MeV compared with the RPA. The SKM* 
result differs from the first two parameterizations. Al- 
ready on the RPA level the strength is less fragmented 
compared to the previous ones. If phonons are included a 
fraction of the narrow RPA peak survives and gives rise to 



FIG. 2: El photo absorption cross section for Sn calculated 
within the present theory with different Skyrme forces. The 
smearing parameter is A = 0.2 MeV. The dotted line are the 
data of ref. [s^. 

a much smaller width than in the previous calculations. 
The method to approximate the strength distribution by 
one single Lorentzian is not well justified in the RPA case 
because the RPA response fragments in two clearly sep- 
arated peaks. If one uses this procedure, one obtains the 
average of the two peaks and a Lorentzian width which 
happens to be large . If the analysis includes the main 
peak only, the RPA width reduces to 1 MeV. In Fig. 

TABLE II: Mean energy and width of the GDR in ^^°Sn and 
^^^Sn calculated within the RPA and the ETFFS. The SLy4 
Skyrme parametrization was used. The El strength is fitted 
to a Lorentzian in the range of (8 - 22) MeV 
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[3] the El strength distributions for io"Sn, ""Sn, ^^°Sii 
and ^'^^Sn are shown. In Table [1111 we compare (Q)RPA 
results with the ones obtained with the present theory. 
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In all cases the energy of the present theory is shifted to 
lower energies by about 2 MeV, independent of the num- 
ber of neutrons, and the widths are increased by 50%. 
As the widths are increased the maximum of the cross 
section is decreased up to 50%. The strength in the two 
closed-shell nuclei ^'^''Sn and ^■^^Sn has about the same 
width, despite the fact that one has added 32 neutrons. 
The two closed-shell nuclei show, as expected, a compact 
and relatively narrow GDR compared to the two open- 
sheU nuclei, "°Sn and ^^osn. 

TABLE III: Integral characteristics of the GDR in "^Sn cal- 
culated within the RPA and with inclusion of phonons with 
three different Skyrme forces. The El strength is fitted to a 
Lorentzian in the (8 - 22) MeV interval. Data are from Ref. 
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We also mention that for ^^^Sn we solved the contin- 
uum RPA. The nearly complete overlap with the dis- 
cretized RPA shows that our discretization procedure, 
together with the smearing parameter A — 200 keV, is 
well justified for the present energy range. As the two 
distributions are hard to distinguish, we do not show the 
continuum result in Fig. [3] 

Our results of ^^°Sn and ^'^■^Sn for the GDR is qualita- 
tively similar to the calculations by Bortignon et al. [s^ . 
They use for their comparison the SIII parametrization 
which gives on the RPA level a centroid energy about 2 
MeV higher than the SLy4 parametrization. Therefore 
SIII has the right tendency. 

However, the far too high incompressibility (see Ta- 
ble IJ) was one of the motivations for the new generation 
of parameterizations with a different density dependence 
{a < 1), which we use in our calculations. Note: there 
is a difference between the definition of the mean en- 
ergy of the total El distribution Eq — mi /mo and the 
mean energy of the GDR which is is used in Table IIIII 
With the first definition the RPA and ETFFS results are 
very similar e.g. for ^^^Sn Eq ^ lAJAMeV (RPA) and 
Eo = \AA7MeV (ETFFS), whereas the mean energy of 
the GDR (Table Hm differ up to 1.8 MeV. If one com- 
pares the theoretical results with the data one has to use 
the latter definition because the experimentalists extract 
their results in the same way from their data. 

Finally we compare in Fig. 2] the original theoretical 
result for the GDR in ^^^Sn and ^^^Sn with the result 
where the isovector force was increased by 20%. Here 
we obtain good agreement with the data. The structures 
beyond 20 MeV are due to non collective ph components 




E [MeV] 

FIG. 3: Strength functions for GDR in "'"Sn, "°Sn, ^^"Sn and 
^^^Sn. The (Q)RPA results (dashed lines) are compared with 
the results of the ETFFS. The SLy4 Skyrme parametrization 
was used. 



which will be smeared out if we include the single particle 
continuum. 



B. Collective electric isoscalar states 

Energetically low-lying collective states, such as the 
3~ mode, are extremely sensitive to the strength of the 
isoscalar particle hole interaction. In Fig. [5] the E2 
strength functions are presented for ^^°Sn and ^■^^Sn. We 
show the results of the RPA and the present theory. The 
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FIG. 4: El photo absorption cross section for ^^"Sn and ^^"Sn 
calculated within the present theory where the isovector force 
has been increased by 20% (full). The dotted lines indicates 
the lorentzian fit to the experimental data [sj, H^. 



FIG. 5: The strength function for 2+ excitation for Sn and 
^^^Sn calculated with (full) and without (dashed) Phonons. 
The smearing parameter is A = 0.2 MeV.The SLy4 Skyrme 
parametrization was used. 



lowest 2+ in ^^°Sn is below zero and, as well, the isoscalar 
giant electric quadrupole resonance is too low compared 
with the data. In ^^^Sn the isoscalar electric quadrupole 
resonance has not yet been measured. In Table IIVI the 
theoretical energies of the low lying collective 2+ and 3~ 
states are compared with the experiments. The RPA re- 
sults for the 2+ modes are in nice agreement with the 
data for both isotopes. The 3~mode in ^'^^Sn is under- 
estimated by 1.35 MeV. In the presence of phonons, the 
collective 3~ modes are pulled down by approximately 
1.5 MeV. The 2+ mode in ^^^Sn is shifted downwards by 
2.36 MeV, and the 2+ mode in ^^"Sn is even shifted be- 
low the ground state energy. Here, however, one has to 
bear in mind that in ^^"Sn the magnitude of the pairing 
gaps plays an important role. 

The reduction of the excitation energies of the collec- 
tive isoscalar modes with respect to the random phase 
approximation has two sources, the modified level den- 
sity in the vicinity of the Fermi level, and the strength of 
the isoscalar particle hole interaction. The compression 
of the level density due to phonons can be understood 
qualitatively by second order perturbation theory: The 
intermediate states generated by the phonon degree of 
freedom have larger excitation energies than the single 
particle level obtained in the mean field approximation. 



Therefore phonons introduce an additional attractive in- 
teraction, whatever the details of the phonon-particle 
coupling are, which increases the level density near the 
Fermi level. The only way to compensate the reduced 
level spacing between particle and hole states is to reduce 
the strength of the isoscalar particle hole interaction in 
comparison with the interaction used in a random phase 
calculation. 

Similar arguments hold for the collective isovector 
modes discussed above. The repulsive isovector parts of 
the particle hole interaction have to be increased com- 
pared with the strength used in the random phase ap- 
proximation in order to compensate the phonon effects. 
This conclusion is independent on details of the theoret- 
ical framework. 



IV. CONCLUSION 

In the present paper we report self-consistent nuclear 
structure calculations performed in the Extended Theory 
of Finite Fermi Systems. Skyrme parameterizations were 
used, which reproduce not only ground state properties of 
nuclei but also excited collective states within the RPA. 
In comparison with RPA calculations, the widths of the 
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giant resonances are brought into better agreement with 
the experimental data. The centroid energies of the res- 
onances are shifted, however. This finding imphes the 
necessity of a renormahzation of Skyrme interactions for 
apphcations beyond the mean field level. Our investiga- 
tion suggests two major renormalizations: a reduction of 
the strength of the isoscalar particle-hole interaction and 
an increase of the isovector strength. 
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TABLE IV: Low-lying collective 2+ and 3" states in ^^"Sn 
and ^^^Sn calculated within and without phonons.The SLy4 
Skyrme parametrization was used. 
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